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Abstract 

The coupling between doubly degenerate electronic states and dou- 
bly degenerate vibrations is analysed for an octahedral system on the 
basis of the introduction of an anharmonic Morse potential for the 
vibronic part. The vibrations are described by anharmonic coherent 
states and their linear coupling with the electronic states is consid- 
ered. The matrix elements of the vibronic interaction are builded and 
the energy levels corresponding to the interaction Hamiltonian are 
derived. 



1 Introduction 

The case of an octahedrally coordinated metal-ion presenting a doubly de- 
generate (non-Kramers) electronic state of symmetry E interacting with a 
doubly degenerate vibrational state of symmetry e constitutes one of the most 



studied Jahn- Teller system [1-9]. In a cluster model the Hamiltonian H 
of the system reads 

H = H, + H, + HjT. (1) 

In Eq. (P, Hf. stands for the Hamiltonian of the electronic part, represents 
the vibrational Hamiltonian and Hjt is the Jahn- Teller interaction. To fix 
the notations, we briefly review the different contributions to H in the usual 
case of harmonic vibrations. 
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The electronic Hamiltonian is well-known. For transition metal-ions, 
it is convenient to work in the strong-field coupling scheme. We denote as 
1^^) and \e) the uncoupled electronic wavefunctions corresponding to ifg and 
transforming according to the two-dimensional irreducible representation E 
(or Eg) of the octahedral (or complete octahedral) point group O (or Oh)- 

We consider vibrations for the doubly degenerate modes with normal co- 
ordinates Qi and Q2 sharing a common frequency u. Let Pi and P2 be 
the generalized momenta associated to Qi and Q2, respectively. Both sets 
{Qi, Q2} and {Pi, P2} transform as basis functions for the irreducible repre- 
sentation E of the group O. The vibrational Hamiltonian is 

H. = ^{P? + P^) + ^{Qi + Ql), (2) 

where m is the effective mass for the two-dimensional harmonic oscillator in 
the plane {Qi, Q2)- The operator can be expressed as 

= htO^Oiai + 02(32 + 1) 

in terms of hcu and of the creation operators aj. and the annihilation operators 
ttk associated to the couples {Qk,Pk) with k = 1,2. 

The coupling between the harmonic vibrations of the ligands and the 
electronic functions localized around the metal-cation can be described by 
means of the linear Jahn- Teller Hamiltonian 

HjT = ^2rmo^EjT{QiDe + Qa^^.). (3) 

In Eq. (^, Dq and are normalized operators in the electronic space with 
the well-known representation [11,12] 

and EjT is the Jahn- Teller energy. The energy Ejt is connected to the 
vibronic coupling constant K by the relationships 

K = ^JhujEjT- 

The interaction Hamiltonian Hjt can thus be rewritten as 
HjT = K[{al + ai)De + (4 + ^2)^^] 
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in terms of K and of the operators a\ and ak with /c = 1,2. A simple 
consideration of the interaction Hamiltonian Hjt shows that the mode Qi 
sphts the two-fold degeneracy of the electronic partners |^) and \e) while the 
mode Q2 mode mixes |^) and \e). 

In order to obtain analytic expressions for quantities of physical interest, 
as for example the energies of the levels or the Ham factor, it is necessary 
to use wavefunctions for the part Hf. + of the Hamiltonian (0). Such 
wavefunctions are products of electronic wavefunctions (of type |^) or \e)) 
with wavefunctions that describe the vibration of the ligands. For the latter 
wavefunctions, we can use coherent states arising from an harmonic oscilla- 
tor with a linear potential that changes the equilibrium point but not the 
frequency. 

The use of the coherent states in the Jahn- Teller effect goes back to 
the idea that the coupling between local vibrations and electronic orbitals 
involves an operator linear in the normal modes in a first approximation 
[13,14]. In this direction, the study of the Jahn- Teller effect using Glauber 
coherent states as well as some applications to particular systems are well- 
known [15-19]. The Glauber coherent states can be written as 



where |00) is the usual ground state for the two-dimensional harmonic oscil- 
lator with vanishing occupation numbers n\ and n2- In Eq. (1), we use the 
vectors p = (pi, P2) and = (a|, 02). In view of 



the Glauber coherent state \G{p)) is an eigenstate of the annihilation oper- 
ator Ofc with the complex eigenvalue pk (A; = 1, 2). 

In orbital doublet systems, the potential energy surface, in the linear cou- 
pling approximation, has a continuous set of minimas which are commonly 
described as forming a "Mexican hat". In actual (i.e., physical) E ® e sys- 
tems, isolated minimas arise either from an ion - lattice interaction term 
which is quadratic in the displacements of the normal coordinates or from 
vibrations which are anharmonic in such coordinates. Both terms warp the 
potential energy surface so that three minimas are formed as a result of one 
of these terms acting either separately or together with the others. For these 





a,|G(p))=p,|G(p)) 
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reasons the surface that represents the adiabatic potential warps and, along 
the bottom of the trough of the hat, three wells occur, alternating regularly 
with three humps known as "tricorns" [8]. Opik and Pryce [2] and O'Brien 
[6] took into account the quadratic terms of vibronic coupling in the Jahn- 
Teller Hamiltonian and established the influence of these terms on the energy 
levels of the system. 

In general, the replacement of the harmonic oscillator potential by the 
anharmonic Morse potential [20,21] turns out to be an efficient way for de- 
scribing anharmonicity of the vibrations. On the other hand, the anharmonic 
coherent states for the Morse potential have been recently investigated with 



great details [^. Therefore, it is the aim of the present paper to consider a 
deviation from the harmonic approximation in the nuclear motions, i.e., to 
explicitly introduce the vibronic anharmonicity by using the Morse potential. 
More precisely, we study in this work the case of an octahedral Jahn- Teller 
system for which the doubly degenerate electronic states are coupled with 
the doubly degenerate vibrational states taken as the anharmonic coherent 
states for a two-dimensional Morse potential. 



2 Morse anharmonic vibrations 

The Hamiltonian H of the E ® e octahedral Jahn- Teller system has the 
same form as in (|l]) but the contribution represents now the vibrational 
Hamiltonian for a two-dimensional Morse oscillator. We thus replace the 
Hamiltonian (^) by the double Morse oscillator Hamiltonian: 

2 

k=l 

with 

2 

^^-fc = |^ + K)(e-'""^-2e-"^^), 

where m is as before the effective mass for the oscillator, a (with o; > 0) 
the anharmonic constant and Vq a positive constant. Furthermore, Xk is the 
displacement from the equilibrium position for the A;-th dimension and pk 
the associated momentum {k = 1,2). (The new conjugated variables Xk and 
Pk correspond to the variables Qk and Pk of Section 1, respectively.) We also 
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use the reduced parameter 



V 



The eigenvalues Eo^Uk) of the operator H^^ are given by 



where 

and 
with 



hQ = A 



rifc = 0,l,---,Ar 



(As usual, [r] denotes the entire part of the real number r.) 

The eigenfunction ip^,. : R+ — >■ R of H^^ corresponding to the level E^ijik] 
is given by |^ 



in terms of the new variable 



(5) 



Vk = i^e 



In Eq. (^), the function F is the confluent hypergeometric function, the 
parameter Sk reads 

u- 1 



Sk 



- rik, 



and the factor 



1 



lT{u - Uk) 



T{v - 2nk) V nk\ 
is a normalization constant. Finally, we shall use the Dirac notation 



|ni,n2) = Vni(2/l)^n2(2/2 
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for the vibrational eigenstates of the two-dimensional Morse oscillator. In 
this notation, the relationships 

{ni,n2\n[,n'2) = 6{n[,ni)S{n2,n2). 

expresses the orthonormalization of the states \ni,n2) on the Hilbert space 
L^(R^, dxidx2). 



3 The anharmonic coherent states 

In some previous works [2^, we have discussed the role played by a creation 
operator an annihilation operator fe_ and an energy operator bo for de- 
scribing the dynamical algebra of the one- dimensional Morse oscillator. The 
results of Ref. can be extended to a two-dimensional isotropic Morse 
oscillator by introducing creation operators 6+^, annihilation operators b_k 
and energy operators bok (with k = 1,2) that parallel the operators 6+, b- 
and bo, respectively. The operators b±k and bok are defined by the differential 
forms 

b±k = i2skTl)^±— -Tt;^ (6) 

dyk Vk 2 



dyi dyk yk 2 2 
with = 1,2. They obey the commutation relations 



[b+k, M = 2bokS{k, I), [b±k, M = ±b±kS{k, I), 

where k and / = 1,2. Therefore, the set {b±k,bok '■ k = 1,2} generates the 
Lie algebra of the group SU(1, 1) ® SU(1, 1). The action of this set on the 
space L'^ili'^, dxidx2) is defined through 

b+i\ni,n2) = \J {rii + - rii - l)\ni + l,n2), 

b+2\nun2) = \l {n2 + l)(z/ - na - l)|ni, + 1), 
b-i\ni,n2) = -\Jniiy - ni)\ni - l,n2), 
b-2\ni,n2) = -yn2(jy^^n2)\ni,n2 - 1), 
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nk)\ni,n2), k = 1,2. 



The operators b±k and 6ofc given by (H) and are s^-dependent, i.e., 
energy-dependent. It is hence necessary to introduce new generators of the 
Lie group SU(1, 1) ® SU(1, 1) which do not depend on Sk and which have an 
action on the states |ni,n2) similar to the one of b±k and bok- This may be 
achieved with the aid of two auxihary variables of a phase type {^k £ [0, Stt] 
for A; = 1, 2). To be precise, let us define the operators 



2_d_ 

i d^k 



oyk Vk 



I d (2 d 



i (9^k V i ^6 



V 



and 



OOfe 



i 

with k = 1,2. The commutators of the operators a+fc, U-k and a^k are 

[a+fc, a-i] = 2aok5{k, I), [a±k, ao/] = ±a±k5{k, I), 
with k and I = 1,2. In addition, we have 

[aoA:,e±^«'] = ±e±^«''(^(A;,/), 

with A; and / = 1,2. 

As a result of our transformation from the 6's to the a's, it appears that 
the wave-function ipn^ belongs to a ray defined via 



^nkiVk.^k 



i^ndVk 



The function is of course an eigenfunction of Hy^ with the eigenvalue 
Eo{nk) and the action of the operators a±k and aok on is 

a+i^nAVki^k) = \J {rik + l){u-nk- l)<l>„,+i(?/fc, ^k)S{k, I), 

a-i^rikiVk^ik) = -\l'nk{y - nk)^nk^i{yk, Ck)S{k, I), 

aoi^nk{yk,Ck) = Sk^nk{yk,Ck)S{k,l), 

where k and I = 1,2. We close this section with a remarkable result: The 
vibrational Hamiltonian can be written as 

Hy = -hil{ali + ttga) 



in terms of the energy operators Cqi and ao2- 
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4 The Jahn- Teller interaction 

The Jahn- Teller interaction Hamiltonian Hjt for the studied system has the 



form 14 



Hjt = + X2), (8) 

where and fi represent the creation and annihilation operators for the 
doubly degenerate electronic states of the system (the label E indicates that 
the operator acts on the electronic states) and k is the strength of the Jahn- 
Teller coupling {k corresponds to K in units of hQ). 

We want to write the operator Hjt in terms of the operators a±k and 
Oofc- For this purpose, we start from the relations 



d 

dyk 



1 00 . 
--E(2ao,r e+^^^ 

^ 771=0 



CL+k 



-ire-*«^a_,, + -fl + 



and 



1 °° r V ^ 

- = -Yl (2aofc)'"-' e+i«'=a+fc + (-l)'"e-'«'=a_fc + -[1 - (-I)'"] 

yk n ^ Z J 



yk m=0 

The coordinate 



Xk = -(lnz/-ln?/fc) 

a 



can be expanded in the form of a series 



Xk = 

« o^^in\ yk, 



or more precisely 



Xk = -\lnu-J2- 



n=l 



n 



2m 



m=0 



+ Yl (2aofc)'"-^[e+'«'=a+fc + (-l)-e-^«'^a_,] 

m=0 

As a result, the Hamiltonian i/j^ can be expressed as 



1 



2 00 



Hjt = }2\niy-YY- 



a 



k=l n=l 



n 



m=0 



2m 
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-k] 



(9) 



m=0 

The latter development in powers of the operators a±k and a^k (with k = 
1,2) is of central importance for the treatment of the vibronic interaction. 

Following the approach in Ref. ||14|, the wavefunctions of the electron- 
phonon system will be written as the overlap 

\/3nz) = r dip\P)e'^^e'^^+ia+ - KrMyu^i)My2.^2), (10) 
Jo 

where /3 refers to the lower branch (/3 = /) or to the upper branch (/3 = u) 
for which we have 

|/) = cos||0)-sin||£), (11) 

|w) = sin||^) + cos||£), (12) 

and where the creation operator a_|_ is defined by 

a+ = a+i cos + a+2 sin (13) 

The phase factor ip in Eqs. (p!0|-p!3|) is an arbitrary angle and the parameter 
z in Eq. ([T0|) characterizes the type of vibronic coupling. For example, the 
case z = ^ corresponds to the states accessible by electric dipole radiation 
from the zero-phonon ground state. 

The states \f3nz) can be rewritten in a new form by expanding the right 
hand-side of (p^O]). This leads to 



p=Oq=0 



EE(-ir— ^C';i|/3,n + p-g,^)c 

p=0 g=0 P- 



= j:i-^r—-^L^n-''{^')\f3pz)o, 

where is a binomial coefficient, \I3pz)q stands for \(3pz) with k = 0, and 
L^""" is a generalized Laguerre polynomial ||23|. Thus, the states \f3nO) cor- 
respond to the case of a vanishing coupling and |/?n|) to the case of a weak 
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coupling. The states (|T^) are eigenstates of the operators b+k (with k = 1,2) 
and of the Hamiltonian (H). Therefore, they are appropriate for obtaining 
the energy due to the Jahn- Teller interaction. 

To find the energy levels of the system with vibronic interactions, we 
calculate the matrix elements of the Hamiltonian Hjt with respect to the 
states (|1^) and use standard perturbation theory. For this purpose, we use 
the notation \n)^ = \(3nz) for the state vectors in the zero-th order approx- 
imation and consider the Jahn- Teller interaction as a perturbation of the 
system described by + Hy. The wavefunctions of the Hamiltonian (|]), in 
the first-order approximation of perturbation theory, have the expression 

n , t,n[m) — hjc)[n) 

As a result, the first-order energy of the system is 

E = Eq + Ejt, 

where 

Eo = ''{n\He + H,\nf 
is the energy without the Jahn- Teller interaction and 

Ejt = °(n|i/jT|n)° 

is the energy of the Jahn- Teller interaction. The calculation of the wavefunc- 



tions and the energy corrections using the coherent states (W) is difficult 
because the expressions for the matrix elements of Hjj- are very compli- 
cated. Moreover, the utilization of the operator exp(ifi;a+) instead of the 
unitary operator exp(i(fi;a+ — K*a)), where a is the adjoint of a+, in building 
the coherent states is rigorous only for the case of the harmonic oscillator 
corresponding to the ordinary Weyl-Heisenberg algebra. For this reason, we 
calculate the matrix elements Ejt by using the Jahn- Teller Hamiltonian (|^). 
It can be proved that the matrix elements corresponding to the last term of 
Eq. (P) vanish. As a final result, we obtain 

1 N \2 p+q-2n 

Ejt = -M Y: i-ir • , , Ll-{^')Ll-{^') 



10 



I 2 oo -1 



k=l 1=0 



i + ^E(^-i-2p) 



2m 



m=0 



where is such that + 1 is the total number of states for the one- 
dimensional Morse oscillator. 



5 Conclusions 

We studied the vibronic coupling between doubly degenerate vibrations and 
doubly degenerate electronic states of an octahedral system. In order to 
describe the vibrations of the system, we used a new type of anharmonic 
coherent states. We described the dynamical system by introducing two 
auxiliary variables (with /c = 1, 2) of extra-phase type and new eigenfunc- 
tions of the Hamiltonian for the anharmonic vibrations. The eigenfunctions 
of the Hamiltonian of the system in absence of the vibronic interaction were 
builded using the overlap of electronic states and these new anharmonic co- 
herent states. In the first-order approximation of perturbation theory, the 
energies due to the Jahn- Teller interaction were calculated. The results show 
that the anharmonic effects on the Jahn- Teller interaction can be expressed 
in a rigorous form and are contained in the constants {k — 1,2). The 
work reported in this paper constitutes the most general algebraic analysis 
of the Jahn- Teller interaction for E®s octahedral systems with anharmonic 
vibrations. 
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